Enhanced ancilla— assisted calibration of a measuring apparatus 
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A quantum measurement can be described by a set of matrices, one for each possible outcome, 
which represents the probability operator-valued measure (POVM) of the sensor. Efficient protocols 
of POVM extraction for arbitrary sensors are required. We present the first experimental POVM 
reconstruction that takes explicit advantage of a quantum resource, i.e. nonclassical correlations with 
an ancillary state. POVM of a photon-number-resolving detector is reconstructed by using strong 
quantum correlations of twin-beams generated by parametric downconversion. Our reconstruction 
method is more statistically robust than POVM reconstruction methods that use classical input 
states. 
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Measurements are at heart of scientific method, be- 
cause they allow to gauge observables in experimental 
tests, leading either to the confirmation or to the ruling 
out of the scientific hypothesis. In quantum mechan- 
ics measurements play a critical role because they con- 
nect the abstract description of quantum phenomena in 
Hilbert space to observable events. In the process of mea- 
surement, a quantum mechanical object interacts with 
a measurement device, and a measurement outcome is 
a result of such interaction. A complete quantum me- 
chanical description of a measurement device is its posi- 
tive operator- valued measure (POVM). In the quantum 
realm, sensor calibration corresponds to determining its 
POVM. In the last decade, the rapid development of in- 
novative quantum technologies promoted POVMs from 
being an abstract theoretical tool to the experimental 
realm. In particular, precise and fully quantum char- 
acterization techniques for sensors [IrH] play a critical 
role for the implementation of quantum information pro- 
cessing, metrology and imaging [STfTtj]. as well as tomog- 
raphy of states [T7H2"I] and operations [25-30 . Quan- 
tum sensor characterization can be thus seen as a si- 
multaneous measurement of multiple parameters, there- 
fore the efficiency of such measurement is of utmost im- 
portance. However, POVM extraction has been experi- 
mentally pursued by brute force methods so far, i.e. by 
probing sensors with a suitably large set of interrelated 
input signals, classical states, yielding slow convergence 
[21 [3]. It was shown [1] that taking advantage of quantum 
resources, e.g. entanglement, can improve convergence 
beyond the traditional methods. Here, we present the 
first experimental POVM's reconstruction that explicitly 
uses a quantum resource, i.e. nonclassical correlations 
with an ancillary state [31] . Our experiment represents a 
major step forward towards quantum mechanical treat- 



ment of sensors: it demonstrates reconstruction of an 
inherently quantum measure of an arbitrary detector's 
performance-its POVM-by realizing for the first time the 
method of ref.[T]. 

A POVM is defined as a set of operators (matrices) n„ 
that give the probability of the measurement outcomes 
via the Born Rule p n = Tr [^n„], where g is the density 
operator describing the system being measured. In prin- 
ciple, it is possible to extract a POVM of a photon detec- 
tor using classical states of light (e.g., coherent states [5]) 
by inverting the Born Rule, after collecting data for a suf- 
ficiently large set of states. A direct inversion, however, 
is a rather delicate and mathematically unstable proce- 
dure, so that even a small uncertainty due to a finite 
statistical sample size can result in a large uncertainty 
in POVM matrix elements. Having a quantum source 
producing on-demand Fock states with a defined photon 
number would simplify the problem significantly, improv- 
ing accuracy in the same measurement time by at least 
of a factor y/N, where N is the number of possible mea- 
surement outcomes for a detector. Unfortunately, there 
are no ideal sources of photon number states. A measure- 
ment scheme based on non-classically correlated bipartite 
systems (beyond TV = 1) is an attractive alternative that 
realizes the full potential of the original scheme of |32j . 
In this case [T], one beam is sent to the Detector Under 
Test (DUT) and the other (an ancilla state) to an ideal 
photon- number-resolving (PNR) detector (i.e., 100% ef- 
ficiency and full photon number resolution), playing the 
role of what in the following we will address as quantum 
tomographer. In this case, by using twin beams, one pro- 
duces heralded (but not pre- defined) Fock states, thus 
yielding a measurement speedup at least y/N. This alter- 
native retains the statistical reliability advantage of the 
on-demand Fock state source. Even with an imperfect 
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tomographer (i.e., efficiency < 1 and no photon-number 
resolution), significant advantage over classical measure- 
ments can be retained. Thus, ancilla-assisted quantum 
schemes, where nonclassical correlations play a key role in 
improving both precision and stability, represent a prac- 
tical advantage of quantum-enabled measurements over 
their classical counterparts. 

Here we provide the first experimental implementation 
of this novel paradigm, and demonstrate an effective re- 
construction method for an arbitrary detector' POVM, 
thus giving the full quantum characterization of its per- 
formance. 

Let us assume that a bipartite system may be prepared 
in a given state, described by the density operator qr, 
and that, besides the measurement made by the detec- 
tor to be calibrated, a known observable with a discrete 
set of outcomes is measured at the tomographer. In our 
scheme, the DUT is a phase-insensitive PNR detector, 
which represents one of the most critical components in 
quantum technology. The detector's POVM elements are 
diagonal operators in the Fock basis, and may be written 
as n n = J2m ^-nm\ m )( m \i where the n„ m 's represent the 
probability of observing n counts when m photons are 
incident on a PNR detector (with the obvious constraint 
that Y) Il nm = 1). n„ m 's are the matrix elements to be 
reconstructed by our measurement. 

In our experiment, the bipartite state consists of 
the optical twin beams qr = \R))((R\, \R)) = 
J2 m -Rml m )l m )i where |m) is the state of one beam with 
m photons, the tomographer is a simple yes/no detec- 
tor with a selectable efficiency 77, defined as including all 
optical losses, and assumes that the detector is live and 
ready to sense incoming light, and R m is the probabil- 
ity amplitude of a particular |m) state. An experimental 
event is a detection of n photons at the DUT paired with 
a measurement outcome ( "yes" or "no" ) at the tomogra- 
pher, which occur with probabilities 

p(n,yes) = ^ n„ m |i? m | 2 [l - (1 - r]) m ], and 

m 

p(n,no) = ^n„ m |i? m | 2 (l-77) m , (1) 

m 

respectively. Upon collecting data to determine p(n, yes) 
and p(n, no), one may invert these relations and recover 
the unknown matrix elements H nm |31j . The distribu- 
tion |i? m | 2 of the bipartite states is determined from the 
photon distribution of the beam addressed to the tomo- 
grapher, which is identical to its twin that is sent to 
the DUT. In this case, the data are the unconditional 
tomographer click events, which occur with probability 
p(no) = J2 m \R m \ 2 (l—r]) m , and allow reliable reconstruc- 
tion of the |i? m | 2 's 22 after collecting data at different 
system detection efficiencies. Note that this procedure 
is much simpler than full quantum tomography |17H20j . 
as no additional calibration is needed to determine the 



|i? m | 2 coefficients, other than the calibration of the effi- 
ciencies at the tomographer. Notice also that entangle- 
ment is not needed to achieve this POVM reconstruction 
of a PNR detector. It is instead, the strong nonclassical 
correlation that enhances the accuracy and stability of 
the reconstruction, thus highlighting the role of squeez- 
ing and ancilla states as a crucial technical resource for 
the development of photonic quantum technologies. 




FIG. 1: (color online) Experimental setup: LiI03 crystal 
pumped with a pulsed 400 nm beam created through second 
harmonic generation (SHG) produces two correlated beams. 
One is sent to the tomographer (T), while its twin is sent 
to the DUT. The tomographer efficiency is varied by rotating 
the linear polarizer. Interference filters (IF) with 20 nm band- 
passes are used to limit out-of-band light on the detectors. An 
FPGA is used for real-time processing and data acquisition. 
The DUT (inset) is a PNR detector made of two Si-SPADs 
connected through 50:50 fibre beam splitter 

The experimental setup (Fig. [T]) consists of an 800 nm 
mode-locked laser, with a repetition rate of 76 MHz, dou- 
bled via second harmonic generation (SHG) to 400 nm, 
which pumps a LHO3 crystal to produce degenerate, but 
non-collinear, photons using parametric down conversion 
(PDC) with Type-I phasematching [5], One of the beams 
from this crystal is sent to the Tomographer, consisting 
of a calcite polarizer (that allows changing the the detec- 
tion efficiency), an interference filter (with a passband of 
20 nm, full width at half maximum) and a silicon Single 
Photon Avalanche Diode (SPAD). The beam is delivered 
to the SPAD through a multimode fiber, which defines 
the spatial collection of the light. Because the down con- 
verted photons have the same polarization in both arms, 
the polarizer can be used to variably attenuate the input 
beam and hence change the efficiency of the tomographer. 
The other PDC beam is directed to our PNR DUT, a 
detector tree consisting of two Si-SPADs, through a cou- 
pling system similar to the tomographer path (i.e., an 
interference filter and a fibre coupler). This two-SPAD 
DUT is able to discriminate between 3 possibilities: 0, 1 
and 2-or-more photo-detections per pulse. With event 0, 
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neither SPAD clicks. With event 1, either SPAD clicks, 
but not both. With event 2, both SPADs click. The 
outputs of the two Si-SPADs of our PNR detector, to- 
gether with the tomographer output and a trigger pulse 
(from the laser), are sent to a Field Programmable Gate 
Array (FPGA)-based processing and data collection sys- 
tem. We distinguish the three possible outcomes of the 
DUT along with the results of the tomographer mea- 
surements. Because detectors have deadtime, the FPGA 
is programmed to avoid taking data when either of the 
detectors in the system is not ready. Before data ac- 
quisition, the tomographer arm polarizer is calibrated to 
provide 20 different system efficiencies needed for the ex- 
periment. 




FIG. 2: (color online) (a) A linearized Poisson distribution 
with respect to detection efficiency. The best fit (red line) 
of the p(no) data (black points) yields a Poisson distribution 
with /i — 0.5983 ± 0.0017 mean photons per pulse, (b) The 
reconstructed bipartite state |-R m | 2 distribution, compared to 
a Poisson distribution with the photon number determined by 
the fit in (a). Uncertainties shown represent the la variations 
in the reconstructions performed on 30 different data-sets. 

To reconstruct the POVM of our DUT, we first deter- 
mine the relative frequencies /(0), /(l) and /(2), respec- 



tively from the number of 0-, I- and 2-click events nor- 
malized to their sum. We also determined the relative 
frequencies of conditional events, paired with tomogra- 
pher's clicks (/(yes|0, r) v ), /(yes|l,77„) /(yes|2, r) v )) and 
no-clicks (/(no|0, rj v ), /(no|l,jj„), /(no|2, r] v )) for each 
efficiency r\ v . As mentioned above, the preliminary step 
in obtaining the POVM elements is the reconstruction 
of the photon number distribution |i? m | 2 22 of the bi- 
partite state. Fig. [2ja) fits the p(uo) data to a Poisson 
distribution with /x = 0.5983 ± 0.0017 mean photons per 
pulse. This is then used to reconstruct the bipartite state 
|i? m | 2 distribution seen in Fig. [2^b) . The experimentally 
reconstructed photon distribution is in excellent agree- 
ment with the Poisson distribution, with a fidelity larger 
than 99.4% (here and in the following, we use the conven- 
tional definition of fidelity as the sum of the square root of 
the product of the experimental and the theoretical prob- 
abilities [22]). Data are shown only up to to = 5 photons 
since in our experiment the probability of observing more 
than 5 photon pairs per pulse is negligible (less than 4 
xlO -4 ). We then substitute the reconstructed |i? m | 2 's 
together with the set of calibrated efficiencies {tj v } into 
Eq. (JIJ , and reconstruct the quantities H nm using a reg- 
ularized least square method [2j |3] to minimize the de- 
viation between the measured and theoretical values of 
the probabilities. In particular, for each output n of the 
DUT, we minimize the deviation between the observed 
Pexp(n, yes) = /(n)/(yes|rz., r] v ) and theoretical probabil- 
ities p(n, yes) if an event n coincided with a click on a 
tomographer and between p exp (n, no) = /(n)/(no|n, r] u ) 
and p(n, no) if an event n was not correlated to a click of 
a tomographer. This is done for each r\ v . 
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FIG. 3: (color online) Reconstruction of the POVM elements for photon numbers up to m = 5. Experimental reconstructed 
(yellow) and theoretical (black) histograms for (a) Ilom, (b) IIi m , and (c) fbm- Quality of reconstruction of POVM elements 
with m < 5 is independently confirmed by observed fidelities of 99.9%. As expected, the accuracy starts deteriorating for input 
states with m > 5. The uncertainty bars represent the statistical fluctuations in the reconstructions performed on 30 different 
data-sets. 

I 

The reconstructed Ho m , IIi m , H2 m are presented in Fig. [3] for input states with up to to = 5 photons. For 
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the first five values (i.e. m < 4), the high fidelities (larger 
than 99.9 %) and low uncertainties highlight the excel- 
lent agreement between theoretical and experimental re- 
sults. The quality of the POVM reconstruction rapidly 
decreases for m > 4, because of the lack of high pho- 
ton number events, as discussed in connection with Fig. 
[H Note that this limitation is not inherent to our cali- 
bration method. In practice, estimating the probabilities 
with sufficient accuracy in the photon number range of 
interest in a finite measurement time requires a bipartite 
state with enough Fock states in that range: our twin 
beam source produces enough states up to m = 4. 

To assess the reliability of the reconstruction, we 
compare the measured probabilities p e xp(^,on) and 
p oxp (n,off) with the ones obtained from Eqs. ([!]) us- 
ing the reconstructed POVM and the reconstructed state 
(see Fig. [I]). The excellent agreement, as seen by the near 
unity fidelities, confirms that the reconstructed POVM 
provides a reliable quantum description of the detection 
process. 
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FIG. 4: (color online) (a) Comparison between measured 
(points) and theoretical (lines) probabilities, p(n, yes), for 
n =1, 2, and 3, for each measurement r\ v . Probabilities for 
n — 2, and 3 are scaled by 5 and 100, respectively. Theoretical 
probabilities are obtained by substituting the measured val- 
ues of the efficiencies r\ v , the reconstructed POVM, and the 
reconstructed |7? m | 2 into Eqs. Panel (b) demonstrates 

the agreement between the theory and experimental data in 
terms of fidelity. 

In conclusion, we have experimentally reconstructed 
the POVM of a photon-number-resolving detector by ex- 
ploiting the quantum correlation of a twin-beam state. 
The reconstructed POVM elements are in excellent agree- 
ment with the theoretically expected ones, as witnessed 
by their fidelities, always above 99.9% for up to four 
incoming photons. Our results represent a major step 
forward towards a quantum photonics for at least two 
reasons. On the one hand, this is the first experimental 
demonstration of an enhanced ancilla-assisted quantum 
detector tomography: we demonstrated the reconstruc- 
tion of a inherently quantum measure of an arbitrary 
detector's performances-its POVM. On the other hand, 
in view of the development of novel PNR detectors with 
improved efficiency, timing jitter, and dynamic range, we 
expect a dramatic growth in the demand of robust, reli- 
able, and fully quantum characterization methods, with 



emphasis on those exploiting quantum resources to go 
beyond the limits of classical measurement. 
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